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Abstract: We show that when the relevant CP violating interactions in leptogenesis are 
fast, the different matter density asymmetries are determined at each instant by a balance 
condition between the amount of asymmetry being created and destroyed. This fact allows 
to understand in a simple way many features of leptogenesis in the strong washout regime. 
In particular, we find some non-trivial effects of flavour changing interactions that conserve 
lepton number, which are specially relevant in models for leptogenesis that rely heavily on 
flavour effects. 
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1. Introduction 

Leptogcncsis is one of the most attractive mechanisms to explain the origin of the matter- 
antimatter asymmetry of the Universe [1]. This is so because it arises naturally in simple 
extensions of the standard model (SM) which can also explain why the neutrino masses are 
so tiny. In this mechanism a lepton asymmetry is produced in the out of equilibrium decay 
of heavy Majorana neutrinos, which is then partially converted into a baryon asymmetry 
by non-perturbative sphaleron processes (see [2] for a complete review). 

Among the three Sakharov conditions [3] for dynamically generating a baryon asym- 
metry, the one requiring out of equilibrium evolution is the least trivial regarding its quan- 
titative aspects. This condition is a consequence of CPT and unitarity, which imply that 
in a universe with fixed temperature and with all the particles having equilibrium phase 
space distributions, there can be no generation of asymmetry (see e.g. [4, 5, 6]). How- 
ever, if there are CP violating interactions mediated by particles which live for some time 
and the temperature of the universe is continuously changing, it is possible to generate a 
matter-antimatter asymmetry. This is so because the thermal bath in the universe is differ- 
ent when the mediating particles decay from what it was when the particles were created. 
Therefore the unitarity relations among CP violating processes, which hold at the level 
of probability amplitudes, cannot be extended to relations among the reaction densities 
-which are the quantities that drive the evolution of the density asymmetries- because they 
also depend on the fluxes of the particles (see [7] for a detailed explanation). In this work 
we will give some insights for the situations in which the CP violating interactions are very 
fast before decoupling and therefore the lifetime of the mediating particles is considerably 
smaller compared to the time scale set by the expansion rate of the universe. 

We will consider the evolution of the different matter density asymmetries as described 
by the classical Boltzmann equations (BE). In the equation for the lepton asymmetry, the 
role of the Yukawa interactions of the heavy neutrinos can be separated in two different 
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type of terms. On one hand, the so caUed "source terms" involve the production of lepton 
asymmetry independently of the amount that is present in the Universe. On the other 
hand, there are "washout terms" proportional to the lepton asymmetry which tend to 
erase the asymmetry populating the Universe at each instant. The main point raised in 
this work is that when the Yukawa interactions are much faster than the expansion rate 
of the Universe, at each instant the thermal bath has time to relax, in the sense that the 
amount of asymmetry being produced equals the amount being destroyed. This observation 
leads to a simple and accurate approximation for the evolution of the lepton asymmetry 
valid whenever the relevant Yukawa interactions are fast relative to the expansion rate 
and the lepton asymmetry is not changing abruptly (we will be more precise about this 
later). The approximation cannot describe the very last stage of leptogenesis when the 
asymmetry freezes, therefore it is not so worthful for obtaining accurate values of the 
final matter-antimatter asymmetry. Notwithstanding it is useful to understand in a simple 
way many features of leptogenesis, in particular the dependence of the produced lepton 
asymmetry on the different parameters of the model. 

The paper is organized as follows. In Sec. 2 we describe the main ideas and intro- 
duce the approximation, starting with the unflavoured case and then generalizing to the 
flavoured one. In Sec. 3 we apply the approximation to study the role of flavour changing 
interactions (FCI), i.e. interactions between different lepton flavours which conserve total 
lepton number, showing some non trivial effects and in Sec. 4 we summarize. 

2. The strong washout balance approximation 

The idea developed here is very general, but for the sake of concreteness we will start 
working in the framework of the type I non-super symmetric model for leptogenesis. In this 
model the particle contents is that of the Standard Model (SM) plus a number of heavy 
Majorana neutrinos Ni. The Ni are SM singlets and only interact with the Higgs and the 
SU(2) lepton doublets through Yukawa interactions. Therefore, in a basis in which Ni are 
mass eigenstates the Lagrangian reads 

C = CsM + iNi^Ni - MiNiNi - Xai V P^iia - X*^~iaPRNji, (2.1) 

where a,i are family indices {a = e,iJ,,r and i = 1,2, 3, ... ), £a are the leptonic SU{2) 
doublets, h = {h^,hP)'^ is the Higgs field [h = iT2h*, with T2 Pauli's second matrix), 
and Pr^l are the chirality projectors. Furthermore, we will mainly focus on hierarchical 
scenarios, for which the Ni mass (Mi) is much smaller than the masses of the other heavy 
neutrinos and we also assume that the lepton asymmetry generated during A^t>i decays is 
not relevant. Hence, from now on the notation is simplified to = A^i and M = Mi. 

Depending on the intensity of the Yukawa interactions, leptogenesis may occur in 
two very different regimes. The weak washout regime is defined by the condition Ftv ^ 
H{T = M), where Tn = M{X^X)ii/8 vr is the decay width of N and H is the Hubble rate. 
In this case the heavy neutrinos decay far out from equilibrium and the erasure of the 
asymmetry produced in the decay epoch is negligible. On the contrary, the strong washout 
regime [Tn ^ H{T = M)] is characterized by small departures from equilibrium and fast 
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processes which tend to erase the lepton asymmetry. Alternatively, the conditions defining 
the weak and strong washout regime can be expressed as rfi -C 10~^ eV and m 3> 10~^ eV, 
respectively, where fh = {X'^ \)iiv^ /M is the so called effective mass and v = 174 GeV is 
the Higgs vacuum expectation value. 

Intuitively one can expect that when the CP violating interactions involving the heavy 
neutrinos are much faster than the expansion rate of the Universe (i.e. in the strong 
washout regime), the lepton asymmetry at each instant takes the value that enforces a 
perfect balance between the production and destruction rates of asymmetry. Actually, the 
requirements for this balance condition to rule the evolution of the lepton asymmetry are 
somewhat more restrictive. It does not hold at the initial stage of leptogenesis, because 
at the beginning the evolution is driven mainly by the source term. Also, when there is 
a change of sign in the lepton asymmetry, its evolution has nothing to do with a balance 
between source and washout. Next we study when and how the balance condition applies 
to the case of unflavoured leptogenesis and then we generalize to the flavoured one. 

2.1 One flavour case 

The BE for the lepton asymmetry has the form 

dYL 



dz 



Siz) - YLWiz) , (2.2) 



where z = M/T, S{z) and YlW{z) are the source and washout terms, respectively, and 
is the lepton number density asymmetry normalized to the entropy density. The balance 
condition introduced above states that, at each instant S{z) YlW{z) and therefore, 
when it is valid, the lepton asymmetry takes the value 

Y.(^) » . (2.3) 

We will call this approximate expression for the lepton density asymmetry the strong 
washout balance approximation (SWBA). 

From the Eq. (2.2) it is clear that the balance condition is equivalent to ^^jj^ ~ and 
that the SWBA will be accurate as long as <C S{z),YlW(z). This helps to understand 
the picture in the regime well described by the SWBA: There are two relevant scales of time, 
a short one given by the rate of the CP violating interactions and a much longer one set 
by the expansion rate of the Universe^ . At each instant there is some production of lepton 
asymmetry (given by the source term) and fast processes which tend to erase it. In the 
short time scale the production and erasure rates are constant and the lepton asymmetry 
takes the value that equilibrates the production and destruction processes. The change in 
the lepton asymmetry is only due to the slow variation of the production and destruction 
rates driven by the temperature decrease of the expanding Universe. It is worth recalling 
that the -slow- expansion of the Universe is fundamental for generating an asymmetry. If 



^Actually, there is yet another important time scale associated with the gauge interactions that keep the 
massless particles (e.g. the leptons and the Higgs) in thermal equilibrium. This is typically much shorter 
than the other scales, as assumed in this work. 
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the Universe were at a fixed temperature, there would be nothing like a source term and 
any initial asymmetry would be erased. Below we comment on an approximate expression 
for the source term that clearly shows its dependence on the expansion rate. 

Considering the validity requirements for the SWBA, namely <^ S{z),YlW{z), it 
is clear why the balance condition does not hold neither at the beginning of leptogenesis, 
nor when the Icpton asymmetry changes sign. In the former situation, Yl is very small and 



dz 



\S{z)\ + 



hence ?a S{z). In the latter, the source term changes sign, therefore 
|y/^iy(z)|. Furthermore, at the end of leptogenesis when the CP violating interactions 
decouple, all the three terms S{z),YlW{z) become very small and the SWBA is 
expected to fail (clearly the evolution during freeze out is not given by a balance condition) . 

In order to see how well the SWBA behaves, let us take first a simple example of 
unfiavoured leptogenesis including only decays and inverse decays. The BE for this case 
are 



ID 



dYN ^ ^_fX^_. 

dz zHs \ Y^^ 

dYB-L I { (Ym a , Yb-l1d 

-^ = ~W-sY\y^^-^)^^^^^ 



(2.4) 



where e parametrizes the amount of CP violation per neutrino decay, Y^ = Ux/s is the 
number density of the quantity x normalized to the entropy density s, Y^-y = Yx — Yy, and 
the superscript eq denotes an equilibrium density. In particular, Yl, Yb, and 1^^' represent 
the lepton density asymmetry, baryon density asymmetry and equilibrium density of one 
relativistic degree of freedom, respectively, all of them normalized to the entropy density. 
Furthermore, 'Jd is the reaction density of decays, which is given by 'Joiz) = ''^n^^'^n, 
with Kn the modified Bessel functions. 

Neglecting the derivative of the Yb-l asymmetry^ leads to the SWBA, 

Yb-l{z) = , (2.5) 

where S{z) = — 1)7d is the source term. 

To use the SWBA as it stands in the above equation, it is necessary to solve the BE 
for the density of the heavy neutrinos because S{z) depends on Yj^j. However, this can be 
avoided using a simple and accurate approximation for S{z). To obtain it, note that the 
source term is equal to e^^j^ (which is also true if scatterings are included consistently [9, 
10]). Therefore, once the density of heavy neutrinos begins to follow closely the equilibrium 



^Here and in the following we work with the B — L density asymmetry, or the B/Z — La asymmetries 

for the flavoured case, because these quantities are conserved by the electroweak sphalcron processes and 
hence the BE are simpler. In any case, all the previous arguments equally apply to these asymmetries. 
Furthermore, when the electroweak sphalerons are fast, the baryon and lepton asymmetries are related to 
Yb-l by Yb — cYb-l and Yl = (c — 1)Yb-l, respectively, with c a constant that depends on the particle 
contents of the model [8]. 
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density, the source term can be approximated by e-gj- and the SWBA reads 

Yb-l{z) = . (2.6) 

zHs 2y/9 

In Fig. 1 we show the exact solution to the BE (2.4) obtained by numerical integration, 
the SWBA using the exact expression for the source term [Eq. (2.5)], and the SWBA with 
the approximated source term [Eq. (2.6)]. For comparison the quantities Y^,S{z), and 
W{z) = 7^ plotted. The effective mass has been set to m = 0.1 eV, which is a 

rather high value, but we have checked that the accuracy of the approximation for effective 
masses one order of magnitude smaller is rather similar (it just works somewhat worse 
at the end of leptogenesis when the Yukawa interactions decouple) . The figure illustrates 
that the SWBA works very well after the lepton asymmetry changes sign and until the 
asymmetry freezes at z ~ 10. 
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Figure 1: Densities and asymmetries in absolute value as a function of z. The thin dashed green 
curve is , the thick dashed green curve is Y/v, the solid red curve is the exact absolute value 
of Yb-l normalized to e, the thick dotted blue curve (labeled SWBAi) is the SWBA for IYb-l/e 
using the exact expression for the source term, and the thin dotted blue curve (labeled SWBA2) 

It is also shown the absolute 



is the SWBA for |FB-i/e| approximating the source term to e- 
value of the source term normalized to e, i.e. (thin dash dotted light blue curve) and the 

inverse decay rate conveniently normalized, —h-^^ = W{z) (thick dash dotted light blue curve). 



We remark that the following things happen at z ~ 10: (I) The washouts decouple, i.e. 
II^(z) becomes less than 1. (II) The density of and the B — L asymmetry (normalized 
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to e) intersect, i.e. |Yb-l| ~ |e|^Ar- The quantity —eYpf has a physical meaning, namely 
it is the total amount of B — L asymmetry that can be produced in the decays of all the 
neutrinos present at a given moment (note that it is equal to the integral of the source 
term from that moment till z oo). Therefore, at z ~ 10 the maximum amount of 
B — L asymmetry that could be produced becomes less than the amount already present. 
(Ill) The lepton asymmetry starts to freeze. This is a consequence of (I) and (II). (IV) The 
SWBA begins to fail. Nevertheless it still works rather fine at this point. To sec this, note 
that when W{z) = 1 the SWBA gives IYb-l/cI = |^^| ^ (the last equality is vahd 
for z > 5), in agreement with (II). 

Prom the figure it can also be seen that the SWBA given in Eqs. (2.5) and (2.6) become 
equal after z ~ 1. This means that the approximation for the source term introduced above 
is very accurate after the population of heavy neutrinos starts to be -almost- in equilibrium. 

To conclude with this first example, we discuss how the SWBA can be used to esti- 
mate the dependence of the matter-antimatter asymmetry on the parameters of the model. 
Besides the CP asymmetry e, the only relevant parameter in this simple example is the 
effective mass rh. Prom the SWBA given in Eq. (2.6) it is clear that Yb-l is inversely pro- 
portional to m during the period ruled by the balance condition. If the temperature 
at which the CP violating interactions decouple were independent of m, then the relation 
Yb-l oc fh~^ would also hold with good accuracy for the final value of the B — L asym- 
metry. Actually T^^c depends on fa, but since the dependence is mild -logarithmic to be 
more precise- and Yb-l decreases only linearly with z, the departure from the "1/m law" 
is not much. This gives an explanation to the known fact that the final B — L asymmetry 
in the strong washout regime depends on fa as l/m^^"*""^, where a has a mild dependence 
with fa (see for instance [11]). Por example, a ~ 0.15 for fa ~ 0.01 eV and it gets smaller 
for larger values of fa. 



2.2 Flavoured case 

When one or more Yukawa interactions of the charged leptons are in equilibrium and have 
higher rates than the Yukawa interactions of the heavy neutrinos, the flavour structure of 
the evolution equations plays an important role [12, 13, 14, 15, 16]. Choosing bases that 
diagonalize the density matrices of leptons and antileptons, the BE for the asymmetries 
can be written as 

^ = S-WYa, (2.7) 

where Ya and S are the column matrices containing the Ya^ density asymmetries (with 
Aq = B/3 — La) and the source terms Sa, respectively (a = 1, . . . , n/ is a flavour index and 
Uf is the number of relevant flavours). Furthermore W is the nj x nj matrix containing 
the washout factors. Note that Sa is proportional to the flavoured CP asymmetry Ca = 

{Ta - ^a)/Y.p(^l3 + f/?), where r„ = r(iV £ah) and r„ = r(A^ IJi) are the decay 
widths of N into Cah and £q,/i, respectively. If all the lepton flavours have fast CP violating 
interactions, the situation is similar to the one described in the unflavoured case, namely 
that at every instant the amount of CP asymmetry produced in each flavour equals its 
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erasure, or equivalently « 0. This leads to the SWBA 

Ya ~ W-^ S, (2.8) 

that generahzes Eq. (2.3) for flavoured leptogenesis. 

Figure (2) shows an example on the use of the SWBA for a flavoured case. This 
time we have included decays, inverse decays and scatterings involving the top quark, but 
for simplicity we have ignored finite temperature corrections for the particle masses and 
couplings [17]. The mass of the heavy neutrino is 10^ GeV and the spectator processes 
active during its decay are also included (see [18, 19] and [15]). The corresponding BE 
are the ones given in [15] but modifying the source term to include the CP asymmetry 
in the scatterings as indicated in [10]. The flavour structure is given by the projectors 
Ka = Tci/^^r^ and by the flavoured CP asymmetries. To give an interesting example 
we have chosen a rather hierarchical set of projectors, namely Kg = 0.03, i^^ = 0.1, and 
K-r = 1 — Kg — = 0.87. Instead, the flavoured CP asymmetries are taken to be all 
equal to some value e, therefore the Aq, asymmetries are simply proportional to e and in 
the figure we plot llAa/el (note that the total CP asymmetry e = Yla^a = 3e). Since 
some of the projectors are very small, the effective mass has been taken considerably large, 
m = 1 eV, in order that the CP violating Yukawa interactions of all the lepton flavours are 
fast (here we are not interested in how much baryon asymmetry is produced but in how 
well the SWBA works in a rather sophisticated example). From the figure it is clear that 
in the strong washout regime, after the change sign and until they begin to freeze, the 
balance condition rules the evolution of the asymmetries also in flavoured cases. 

0.01 
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Z 

Figure 2: An example on the use of the SWBA for flavoured leptogenesis. The solid curves give 
the asymmetries \Ya^ /e| as a function of z obtained by solving numerically the BE, while the dotted 
ones correspond to the SWBA represented symbolically in Eq. (2.8). Also plotted are Yn (thick 
black dashed curve) and Y^^ (thin black dashed curve). 
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3. Application to flavour changing interactions 



In this section we will apply the SWBA to study the role of FCI. This type of interactions 
is important, for example, when the Yukawa couplings of the heavy neutrinos are large, 
so that the processes mediated by virtual heavy neutrinos are fast, and also, as remarked 
in [20], in supersymmetric scenarios due to off diagonal soft breaking slepton masses. They 
can potentially enforce the equality of the lepton flavour density asymmetries, inhibiting 
large enhancements of the final baryon asymmetry due to flavour effects. Their effect could 
be even more drastic in models that have a null total CP asymmetry and rely completely 
on flavour effects (see e.g. [21, 22, 23]). This class of models cannot produce any baryon 
asymmetry if all the lepton flavour asymmetries are equal. Next we will argue that having 
fast FCI is not a sufficient condition for the flavour asymmetries to become equal and, in 
particular, that their effect in the models with null total CP asymmetry is less harmful 
than what could seem a priori. 

Wc introduce the notation 7^ for the sum of the interaction rates that produce |AL| = 
|ALq,| = 1 and 7q,^ for the sum of the ones that conserve L but violate and Lp by one 
unit, all of them normalized in such a way that the BE can be written simply as^ 



dz 
d^A 



dz 
dlA. 



= S^- T^^A, - let. {Ya, - l^Ae) - l^r {Ya, - ^A.) , (3-1) 
=Sr- 7rl^A. - ler (^^A. - ^Ae) - T/xr (^A. - IaJ • 

When all the lepton flavours have fast CP violating interactions, the values of the Y/\^ 
asymmetries are determined by the condition that all the derivatives be negligible. 

The solution to the linear system of equations that results after setting the derivatives of 
the Ya„ asymmetries to zero is 

^Ae = [itilrSe + lenlASe + S^) + lerltiiSe + Sr) + l^ril^ + lT)Se 
+ {lenll.r+leTlnT + lenleT){Se + S,^ + Sr)] / \W\ , 



^A;, = [irleSn + li^rleiSn + Sr) + leiil^S^ + Sg) + lerilr + le)S^ 
+ {l^i.TleT + leiileT + ltiTleii){S^j, + Sr + Se)] / \W\ , 

^A^ = [lelixSr + lerlniSr + Se) + lurleiSr + -5^) + lenile + 1h)St 
+ {leTlen + lnTlen + leTlnT){ST + Se + Si^)] / \W\ , 

where is the determinant of the washout matrix, 

\W\ =lelfilT + It^rleilt, + Ir) + lerlfi{le + Ir) + lenlAle + 7^) 
+ {letiler + lefilfir + leTltiT){le + + Ir) ■ 



(3.2) 



(3.3) 



^Spectator processes can easily be included in the analysis, but for simplicity we will neglect them. 
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Each Aq asymmetry in the above equations is given by a sum of terms which contain 
products of two interaction rates and one source factor, all of them divided by the common 
factor \W\. Some terms involve the product of two AL = 1 interaction rates, others have 
the product of one AL = 1 interaction rate with one FCI rate, and finally there are terms 
with products of two FCI rates. The sum of this last type of terms is the same for every 
Aq, asymmetry, therefore all the Yaq (ct = e, /x, r) will be equal in the case that these terms 
dominate over the others. But for this to happen it is not enough to have fast FCI rates (i.e. 
faster than the expansion rate). Basically two more conditions must be satisfied. In order 
to state them in a simple way, let us assume that all the FCI rates have the same order 
of magnitude, with their average value denoted by 7fci(-2^)- In addition, all the AL = 1 
interaction rates are also taken to be similar in magnitude, with an average value equal to 
7al=i{z). Then, the aforementioned conditions are 

(I) 7fci(^) » 1al=i{z) , 

(II) EaSa^l. 

In particular, the condition (II) for ensuring equality of all the lepton flavour asymmetries is 
not trivial a priori, but follows rather directly from the analysis based on the SWBA. Since 
the models with null total CP asymmetry verify Sa = 0, the asymmetries Y^^ will not 
be in general equal in these cases, even under the presence of very fast FCI (this holds 
during the period ruled by the balance condition but not after the AL = 1 interactions 
decouple, as will be further discussed in the example given below). 

From the Eqs. (3.2) it is also possible to estimate the dependence of the density asym- 
metries on the intensity of the FCI. We take 7pci oc m^^j, where m^^ parametrizes the 
intensity of the FCI (the reason for this notation will become clear below). For the case 
e the Eqs. (3.2) indicate that if 7fci ^ 7al=i before the CP violating interactions 
decouple, then |lAa| oc l/rhg^. Note that this statement, as well as the conditions (I) and 
(II) for having equal flavour asymmetries, hold for generic flavour structures of the Yukawa 
couplings, while special cases like the presence of null projectors or some projectors being 
equal have to be treated separately. 

To show a detailed example we consider the case of soft leptogenesis [24, 25], taking 
into account the full flavour structure [26] and including the FCI induced by off diagonal 
soft supersymmetry breaking masses for the slepton doublets [20]. A difference between 
the example given here and common analysis in the literature is that, in order to get a null 
total CP asymmetry, we do not assume the trilinear couplings of the soft breaking 
terms AaiXaiNi^h to be flavour independent. Hence the flavoured CP asymmetries in the 
resonant condition can be written as ea{T) = KaiToa[Aai] Abf{T)/M (for one generation 
of heavy neutrinos) , where A^^- (T) is a thermal factor that goes to when T — >■ [25] . 
For the example we have chosen Kg = 0.1, = 0.2, and Kt = 0.7, whereas Im[^Q,i] = 
|>l|sin^Q;, with \A\ = 10^ GeV, (f>e = 7r/2, = 7r/4, and (j)r equal to the value that 
makes e = = 0. Furthermore, we have calculated and included the rates of the FCI 

induced by the off diagonal soft breaking slepton masses, which are taken to be all equal to 
a parameter called rhod (see [20] for a detailed description), so that 7fci oc rfi^^. The FCI 
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rates also depend on tan /3, taken here to be equal to 30. For more details on the evolution 
equations and parameters involved the reader is referred to [26]. 

Figures (3) and (4) contain the results of our analysis. In Figure (3) the final baryon 
asymmetry is plotted as a function of m^^ for a fixed value of the effective mass m = 0.2 eV 
and Ml = 10^ GeV. It is clear that the behavior described above is verified. For small 
rhod the asymmetry is independent of rfiod, but when the FCI become fast enough, the 
final asymmetry satisfies Yb-l oc The -approximate- value of the cosmological 

asymmetry is indicated with a horizontal line and it can be seen that in this example the 
cosmological asymmetry is achieved under FCI dominating over all the other rates, even 
though the total CP asymmetry is null. 
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Figure 3: The final baryon asymmetry as a function of m^^ for the example described in the text. 
The horizontal line corresponds to the cosmological asymmetry and the small square is a reference 
for Fig. (4). 

Figure (4a) shows the evolution of the density asymmetries corresponding to the point 
signaled in Fig. (3), which is representative of a situation with FCI much faster than 
the expansion rate and the AL = 1 interactions during most of the decay epoch [see 
Fig. (4b)]. As can be seen the SWBA reproduces with reasonable accuracy the exact 
results. In particular the "predicted" non equality of all the flavour asymmetries during 
the leptogenesis era, even with very fast FCI, is verified. Note that the flavour asymmetries 
do equalize at the end of leptogenesis, but this happens when the AL = 1 interactions are 
decoupling and the B—L asymmetry is freezing, so that Yb-l is not affected much. Instead, 
if the equality among the asymmetries had hold at earlier times, the total asymmetry would 
have been null, since e = 0. 

From Figure (4a) it can also be seen that in this case the SWBA ceases to be good 
somewhat before the freeze out, due to the late change of sign of ^At- Nevertheless, the 
dependence of Yb~l on a given parameter of the model does not change much between the 
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period ruled by the balance condition and the end of leptogenesis, as long as T^^gc does not 
depend on this parameter (remember the discussion at the end of Sec. 2.1 and take T^ec for 
the flavoured case equal to the temperature at which the largest CP violating interaction 
decouples). This is why the SWBA remains useful for obtaining the dependence of the final 
asymmetry on some parameters of the model, in particular on rhod which does not change 
Tdec- On the contrary, the study of the dependence of Yb-l on rh in this example is more 
involved. While the SWBA holds, Yb-l is independent of rh when 7fci ^ 7al=i- But the 
final B — L asymmetry gets some dependence on rh because T^ec decreases with increasing 
rh. Contrary to the simple case described in Sec. 2.1, here Yb-l is decreasing exponentially 
before freezing out^, hence a small change in T^ec produces a big change in the final value of 
Yb-l- It can be shown that if the CP asymmetries were independent of the temperature, 
then the final asymmetry would approximately satisfy Yb-l ^ 1/rh, as a result of the 
combination of the logarithmic dependence of T^ec on rh with the exponential decrease 
of Yb-l before freeze out. But actually the CP asymmetries decrease with temperature, 
which adds an additional dependence of the final B — L asymmetry on m. This could 
be studied more precisely using the SWBA together with a knowledge of T^ec = Toeci'rh), 
however it is not the goal of this work to give more details on this point. 

As a final remark, we want to refer to a model of leptogenesis with null total CP 
asymmetry presented in [22]. Actually, a brief mention of the basic idea developed in this 
work was made there to gain understanding about the dependence of the baryon asymmetry 
on the different parameters. In this model, as well as other TeV scale models for leptogenesis 
[see e.g. [27]], the electroweak sphalerons decouple before the end of leptogenesis. Then, 
the baryon asymmetry often freezes in the regime well described by the balance condition 
and in these cases the results obtained from the SWBA are very accurate. 



*This can also be understood from the SWBA. The Eqs. (3.2) for the case e = and 7fci ^ 7Ai=i 
roughly imply that Ya„,{z) ~ S{z) /^fci{z), where S{z) is some combination of the source terms. Since 

dY"'' _ 

S{z) oc -jf- oc e ^ for a > 1 and 7pci is not exponentially suppressed, the exponential decrease of Yb-l 
with increasing z follows. 
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Figure 4: 

The asymmetries and rates as a function of z, corresponding to tlie point signaled with a small 
square in Fig. (3). (a) The asymmetries cYAei cYa^, and Yb = Yla c^A„ as a function of z, where 
c = 8/23 is the conversion factor between the B — L and the B asymmetries due to the electroweak 
sphaleron processes. The solid curves give the asymmetries obtained by solving numerically the 
BE, while the dotted ones correspond to the SWBA. The asymmetry Ya^ is approximately equal 
to and hence it has not been plotted, (b) The reaction densities of decays summed over all 
flavours (7_d), CP violating scatterings involving the top quark also summed over all flavours (7^), 
and FCI induced by off diagonal soft breaking slepton masses averaged over all flavours (josm), 
as a function of z. All the reaction densities have been normalized to Hni, where ni denotes the 
equilibrium density of one relativistic degree of freedom. 
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4. Conclusions 

It has been shown that when the CP violating interactions in leptogenesis are fast, the 
classical evolution of the lepton asymmetry is ruled by an equality condition between the 
amount of asymmetry that is being generated and erased. This fact can be used to obtain 
a simple approximation for the asymmetry, called here SWBA, which can generically be 
written as Yj, = S{z)/W{z). Here S{z) is the source term and W{z) is a combination of 
interaction rates such that W{z)Yl is the whole washout term. This relation is also valid 
when different lepton flavours participate in the dynamics. In many models the source term 
can be obtained with high accuracy just from the rate of change of the heavy neutrinos 
equilibrium density induced by the expansion of the Universe. Therefore, the SWBA often 
does not require any integration of differential equations. 

To be more precise, the balance condition between the generation and destruction of 
asymmetry holds as long as the CP violating interactions remain fast with the exception 
of the initial period in which the lepton asymmetry is starting to be generated or when 
it is changing sign. More importantly, the SWBA fails at the end of leptogenesis because 
the Yukawa interactions of the heavy neutrinos become slow, therefore in general it is not 
useful for obtaining accurate values of the final asymmetry. However, as explained in this 
work, the SWBA allows to understand in a simple and rather accurate way the dependence 
of the asymmetry on the parameters of the model. In particular, we have shown that it 
is useful for studying the models with null total CP asymmetry and fast FCI. In these 
cases flavour effects are mandatory to have successful leptogenesis and hence FCI could 
be potentially harmful. However, we have found that fast FCI do not necessarily spoil 
baryogcncsis, since the final asymmetry has a quite mild dependence on the intensity of 
these rates, namely Yb oc 1/?77,^^ (with the rate of the FCI being proportional to m^^). This 
relation is comparable to the dependence of the baryon asymmetry on the effective mass 
in simple models for leptogenesis which do not require flavour effects to be successful. 

Acknowledgments 

We thank E. Nardi and J. Salvado for useful discussions and comments. We are specially 
indebted to M. C. Gonzalez-Garcia and N. Rius for discussions and careful reading of the 
manuscript. 

This work is supported by Spanish MICCIN under grants 2007-66665-C02-01, ACI2009- 
1038, and Consolidcr-Ingenio 2010 CUP (CSD2008-00037), by CUR Generalitat de Catalunya 
under project 2009SGR502, and by USA-NSF grant PHY-0653342. 

References 

[1] M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification, Phys. Lett. B174 

(1986) 45. 

[2] S. Davidson, E. Nardi, and Y. Nir, Leptogenesis, Phys. Rept. 466 (2008) 105-177, 
[hep-ph/0802.2962]. 



- 13 - 



[3] A. D. Sakharov, Violation of CP Invariance, c Asymmetry, and Baryon Asymmetry of the 
Universe, Pisma Zh. Eksp. Teor. Fiz. 5 (1967) 32 35. 

S. Weinberg, C'osmological Production of Baryons, Phys. Rev. Lett. 42 (1979) 850 853. 

E. W. Kolb and S. Wolfram, Baryon Number Generation in the Early Universe, Nucl. Phys. 
B172 (1980) 224. 

A. D. Dolgov and Y. B. Zeldovich, Cosmology and Elementary Particles, Rev. Mod. Phys. 53 
(1981) 1-41. 

E. Roulet, L. Covi, and F. Vissani, On the CP asymmetries in Majorana neutrino decays, 
Phys. Lett. B424 (1998) 101-105, [hep-ph/9712468]. 

J. A. Harvey and M. S. Turner, Cosmological baryon and lepton number in the presence of 
electroweak fermion number violation, Phys. Rev. D42 (1990) 3344-3349. 

A. Abada, S. Davidson, A. Ibarra, F. Josse-Michaux, M. Losada, and A. Riotto, Flavour 
matters in leptogenesis, JHEP 09 (2006) 010, [hep-ph/0605281]. 

E. Nardi, J. Racker, and E. Roulct, CP violation in scatterings, three body processes and the 
Boltzmann equations for leptogenesis, JHEP 09 (2007) 090, [hep-ph/0707.0378]. 

W. Buchmiiller, P. Di Bari, and M. Pliimacher, Leptogenesis for pedestrians, Ann. Phys. 315 
(2005) 305, [hep-ph/0401240]. 

R. Barbieri, P. Creminelli, A. Strumia, and N. Tetradis, Baryogenesis through leptogenesis, 
Nucl. Phys. B575 (2000) 61-77, [hep-ph/9911315]. 

T. Endoh, T. Morozumi, and Z. Xiong, Primordial lepton family asymmetries in seesaw 
model. Prog. Theor. Phys. Ill (2004) 123-149, [hep-ph/0308276] . 

A. Abada, S. Davidson, F. X. Josse-Michaux, M. Losada, and A. Riotto, Flavour issues in 
leptogenesis, JCAP 0604 (2006) 004, [hep-ph/0601083]. 

E. Nardi, Y. Nir, E. Roulct, and J. Racker, The importance of flavor in leptogenesis, JHEP 
01 (2006) 164, [liep-pli/0601084]. 

S. Blanchct, P. Di Bari, and G. G. Raffelt, Quantum Zeno effect and the impact of flavor in 
leptogenesis, JCAP 070S (2007) 012, [hep-ph/0611337]. 

G. F. Giudicc, A. Notari, M. Raidal, A. Riotto, and A. Strumia, Towards a complete theory of 
thermal leptogenesis in the SM and MSSM, Nucl. Phys. B685 (2004) 89, [hep-ph/0310123]. 

W. Buchmiiller and M. Pliimacher, Spectator processes and baryogenesis, Phys. Lett. B511 
(2001) 74, [hep-ph/0104189]. 

E. Nardi, Y. Nir, J. Racker, and E. Roulet, On Higgs and sphaleron effects during the 
leptogenesis era, JHEP 01 (2006) 068, [hep-ph/0512052]. 

D. Aristizabal Sierra, M. Losada, and E. Nardi, Lepton Flavor Equilibration and 
Leptogenesis, JCAP 0912 (2009) 015, [hep-ph/0905.0662]. 

D. Aristizabal Sierra, L. A. Munoz, and E. Nardi, Purely Flavored Leptogenesis, Phys. Rev. 
D80 (2009) 016007, [liep-pli/0904.3043]. 

M. C. Gonzalez- Garcia, J. Racker, and N. Rius, Leptogenesis without violation of B-L, JHEP 
11 (2009) 079, [hep-ph/0909.3518]. 



- 14 - 



[23] S. Antusch, S. Blanchet, M. Blcnnow, and E. Fernandez-Martinez, Non-unitary Leptonic 
Mixing and Leptogenesis, JHEP 01 (2010) 017, [hep-ph/0910 . 5957]. 

[24] Y. Grossman, T. Kashti, Y. Nir, and E. Roulet, Leptogenesis from Supersymmetry Breaking, 
Phys. Rev. Lett. 91 (2003) 251801, [hep-ph/030708l]. 

[25] G. D'Ambrosio, G. F. Giudice, and M. Raidal, Soft leptogenesis, Phys. Lett. B575 (2003) 
75-84, [hep-ph/0308031]. 

[26] C. S. Fong and M. C. Gonzalez-Garcia, Flavoured Soft Leptogenesis, JHEP 06 (2008) 076, 
[hep-ph/0804.4471]. 

[27] A. Pilaftsis, Electroweak Resonant Leptogenesis in the Singlet Majoron Model, Phys. Rev. 
D78 (2008) 013008, [hep-ph/0805 . 1677]. 



-15- 



